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Abstract 

We give an algorithm to construct a translation-invariant transport kernel between 
ergodic stationary random measures 4> and T on R d , given that they have equal 
intensities. As a result, this yields a construction of a shift-coupling of an ergodic 
stationary random measure and its Palm version. This algorithm constructs the 
transport kernel in a deterministic manner given realizations ip and t/> of the measures. 
The (non-constructive) existence of such a transport kernel was proved in f8j. Our 
algorithm is a generalization of the work of ED, in which a construction is provided for 
the Lebesgue measure and an ergodic simple point process. In the general case, we 
limit ourselves to what we call constrained densities and transport kernels. We give a 
definition of stability of constrained densities and introduce our construction algorithm 
inspired by the Gale-Shapley stable marriage algorithm. For stable constrained 
densities, we study existence, uniqueness, monotonicity w.r.t. the measures and 
boundedness. 
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1 Introduction 


For a random measure tk on E, d , there are a number of equivalent definitions for 
the Palm distribution of \k (see Section [2J. Heuristically speaking, in the case that >k is 
stationary and ergodic, the Palm version of \k is obtained by viewing \k from a typical 
point of \k. The result of 111311 shows that there exists a random point Y such that by 
viewing tk from Y; i.e. by translating ')/ by vector —Y, we get exactly the Palm version of 
'k. In other words, there exists a coupling of \k and its Palm version such that almost 
surely each one is a translated version of the other one. Such a coupling is called a 
shift-coupling. 

To obtain a shift-coupling, one can use a (random) balancing transport kernel T that 
transports a multiple of the Lebesgue measure to \k, where by a (random) transport 
kernel we mean a (random) function that assigns to each point s £ !R d and w in the 
probability space, a probability measure T u (s, •) on R d . This measure can be interpreted 
as how the infinitesimal mass at s is distributed in the space. Given that T depends on >k 
in a translation-invariant manner (which is called flow-adapted here), then choosing Y 
with distribution T w (0, •) gives a shift-coupling of 'k and its Palm version (see (8) and |5j 
as explained after Theorem 4.91. A case of special interest is when T is a balancing 
allocation in which the measure T w (s, •) is a Dirac measure a.e. In this case, given ik, 
the vector Y defined above is deterministic and the shift-coupling does not need extra 
randomness. 

Let $ and \k be jointly stationary and ergodic random measures on R d . Considering 
flow-adapted transport kernels balancing $ and T has been of great interest recently. 
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The abstract existence of such transports was proved in J8) (which is merely based 
on |13l ) provided that the intensities are equal, finite and positive. Nevertheless, there 
has been tremendous interest in the construction of such transports in special cases in 
the recent years. One of the reasons for this interest is that such transport kernels lead 
to explicit shift-coupling of the Palm distributions of $ and T, as explained above. The 
constructions were motivated by Liggett lllOll . who constructed a balancing allocation for 
an ergodic simple point process in dimension one. The landmark in this topic is (3j which 
generalizes 11011 to arbitrary dimensions. There are several other constructions in the 
literature for the case of simple point processes, such as gravitational allocation m, 
optimal transport ©, one-sided stable allocation on the line O, etc. 

In this paper we give an algorithm that works in the general case and enables us to 

• construct a flow-adapted transport kernel balancing two arbitrary jointly stationary 
and ergodic random measures on R d with equal intensities (Theorem |4.8| ) 

• and constructs a shift-coupling for an arbitrary stationary ergodic random measure 
on R d and its Palm version (Theorem |4.9| l. 

To do this, we generalize the notion of stable allocations introduced in O to what 
we call stable constrained densities, where the notion of constrained densities is a 
special case of capacity constrained transport kernels introduced in 13. The first, and in 
fact deterministic, result (Theorem |4.13D is that stable constrained densities exist and 
one can be given by our algorithm (Algorithm |4.4| l which is inspired by the continuum 
version of the Gale-Shapley algorithm in (3j. Another important result is considering the 
algorithm in the random case described above. Other results deal with monotonicity 
(Theorem |4.2ip and optimality (Corollary |4.22[ l properties of stable constrained densities, 
uniqueness (Theorem |4.24[ ) and boundedness of the support of the mass transported 
from and to a region (Theorem |4.25[ ). These results are in the spirit of the seminal 
papers Q and lO and generalize some of their results. We also introduce the notion 
of Voronoi transport kernel with respect to a measure, which generalizes the notion of 
Voronoi diagram for a discrete set. It helps us in proving some statements, but it can be 
interesting in its own. The construction and results in this paper can be generalized to 
random measures on a locally compact Abelian group and also to non-ergodic cases (only 
equality of sample intensities is important) in the setting of (8j. We don't go through 
these general cases to stay focused on the main ideas. 

We show that using transport kernels is inevitable in the general case by providing 
examples where no flow-adapted balancing allocation exists. However, some results and 
open problems on the existence and construction of flow-adapted balancing allocations 
are addressed in Subsection l4.6l 

The paper is structured as follows. Preliminaries about random measures and 
transport kernels are reviewed in Section [2] Since we generalize some works of (3), in 
Section [3] we review the idea of the algorithm in ||3j as a motivation of our work. Our 
main definitions and results are presented in Section [4] and the proofs are postponed 
to Section [5] In Subsection |4 .1 1 we define constrained transports. We introduce our 
algorithm in Subsection |4.2 1 and state the shift-coupling and balancing properties in the 
random case. In Subsection |4.3| we generalize the notion of stability which is a key tool for 
proving the results. Other properties of stable transports are provided in Subsection |4.4| 
Voronoi transport kernels are defined in Subsection |4.5| In Subsection |4.6| we study the 
existence and constructions of balancing allocations. Finally, some examples are given 
in Section[6]which are addressed in the text. 
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2 Preliminaries 

Let Q be the Borel u-field on R d and be the Lebesgue measure on (R rf ,t/). We 
denote by M the set of all non-negative locally finite measures on (R d , Q), and by M. the 
smallest cr-field on M such that the mappings p >->• p(B) are measurable for all B € Q. 
All measures in this paper are assumed to be members of M. 

Recall that the complement, the interior and the boundary of a set A are denoted 
by A c , A° and dA respectively. We denote The identity function on A by 1^. Also, the 
closed ball with center x € R d and radius r is denoted by B r ( x). 

For a measure ip and a non-negative measurable function / : R d —► R, we denote by 
ftp the measure B f B f(s)<p(ds). In this article, by tpi > tp 2 we mean ipi(B) > ipi{B) 
for all Beg. A weighted transport kernel is a measurable map T : R d — > M. For 
simplicity of notation, we denote (T(s)) (B) by T(s , B). Intuitively, we can think of T(s, B) 
as the portion of the infinitesimal mass at s transported to the set B. Given measures ip 
and ip, T is called non-weighted or Markovian if 

T(s, R d ) = 1 for y>-almost all s £ R d 
and it is called (pp, pi) -balancing if it transports <p to ip; i.e. 


/R d 


T{s, -)ip(ds) = ip{-). 


In this paper, by a transport kernel we mean a non-weighted transport kernel. Being 
non-weighted is equivalent to the condition that the total mass transported from a set B, 
which is f B T(s , R d )</j(ds), is equal to p(B). Being (tp, VO-balancing means that the mass 
transported into B is equal to ip(B). If T is (ip, ip) — balancing transport kernel, then T is 
a Markovian kernel which transports ip to ip. 

A (ip, ip) -balancing allocation is a transport kernel T such that for p>- a.e. s, T(s, ■) 
is a Dirac measure <5 r ( s ) and V 3 ( r_1 (')) = ^(‘)- See Definition 4.31 for a precise definition 


of an allocation. For more details on transport kernels, see [|8J. 

In this work, we fix a measurable space (fl, J 7 ) equipped with a measurable flow 
0 S : fl —>• fl for s € R d ; i.e. (w, s) 8 s uj is measurable, 0 O is the identity on f 2 and 


9 s o e t =9 s+u Vs,tG R d . 


With an abuse of notation, we use 9 S also for natural flows on the space of functions on 
R d or R d x R d , on M, and on the space of weighted transport kernels; i.e. 


O s f{x) 

o s f(x,y) 

9 s tp(B) 

O s T{x,B) 


f(x + s), 
f(x + s,y + s), 
p(B + s), 

T(x + s,B + s). 


Using these conventions, a measurable function F from U to a flow-equipped space is 
called flow-adapted if F(9 s ijj) = 9 s F(co). We denote F(u) by F^ too. 

A random measure is a pair (P, 4>), where P is a probability measure on Q and 
f : O -> M is a measurable function. The distribution of (P, <1>) is the push-forward 
measure P($ _1 (-)) on M. 

A probability measure P on fl is stationary if it is invariant under 9 S for all s in R fi . 
A stationary random measure is a random measure (P, 4>) such that P is a stationary 
probability measure and <1> is flow-adapted. This implies stationarity in the usual sense, 
which is translation-invariance of its distribution; i.e. P[$e A] = P [0 K T € A] for any 
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A £ A4 and s £ R d . Moreover, if (P, 4/) is another stationary random measure (with the 
same P), then $ and T are jointly stationary in the usual sense; i.e. the distribution of 
the pair (4>, 'h) is invariant by translations. When there is no ambiguity about P, we may 
say that $ is a stationary random measure in the sense given above. Therefore, when 
we mention two stationary random measures without definig the probability measure(s), 
they are assumed to be jointly stationary. 

A stationary probability measure P on SI is ergodic if for any event A £ T that is 
invariant under all 9 S , we have P(A) £ {0,1}. The stationary random measure (P, 4>) is 
ergodic if P is ergodic. It is clear that if 4/ is a measure derived from $ in a translation- 
invariant manner, then stationarity (resp. ergodicity) of (P, 4>) implies stationarity (resp. 
ergodicity) of (P, T). For simplicity, we use the term ergodic random measure instead of 
ergodic stationary random measure. 

Since we deal with random measures in this work, we have two different notions of 
'almost everywhere'; one corresponding to probability and the other corresponding to 
the measures on R d . To avoid confusion, we preserve the phrase almost sure(ly) (or a.s.) 
for probability measures on 11 and we use almost everywhere and almost all (denoted by 
a.e. and a.a.) for measures on R d . 

Lemma 2.1. Let V C R d be a subset that contains arbitrarily large balls. If $ is a 
stationary random measure such that a.s. 4>(R d ) > 0, then a.s. 4>(E) = oo. 

The intensity of a stationary random measure (P, 4>) is the unique constant A £ [0, oo] 
such that E [4?(R)] = A Cd{B) for an arbitrary Borel set B £ Q. If A is positive and finite, 
then the Palm version of (P, 4>) is the random measure (P$, 4>) in which 


R$(AL) 


A C d (B) E 



^A{S s )L B (s)^(ds) , 


( 2 . 1 ) 


for all A £ T, where 9 S is interpreted as the random element 9 S (w) and B is an arbitrary 
Borel set B £ Q with positive and finite Lebesgue measure. Note that $ stays the same 
and only the probability measure is changed. It can be shown that stationarity implies 
that this definition is independent of the choice of B and formalizes the intuition of the 
Palm distribution given in the introduction. The above equation is equivalent to the fact 


that 


E$ [H] 


A C d (B) E 



H{9 s )t B (s)<t>(ds) 


( 2 . 2 ) 


for all non-negative measurable functions H : LI —>• R, where E$ is expectation with 
respect to P$. By this equation, we can think of E$ [H] as averaging H(9 s lj) over the 
points s £ B and lj £ LI. 

The refined Campbell theorem states that for all non-negative measurable functions 
H : LI x R d R, 


AE<;> 


H(9q, s)ds 


= E 


H(0 s ,s)$(ds) 


(2.3) 


where ds is a short form of Cd(ds). If (P, \&) is also a stationary random measure on the 
same space, we have Neveu's exchange formula 


A$E$ 


H(9 0 ,s)^>{ds) 


= A\pEg> 


H(9 S , —s)$(ds) 


(2.4) 


where A$ and A$ are the intensities of 4> and 4' respectively. Interested readers may refer 
to 11121 sections 3.3 and 3.4, for more details on the properties of Palm distributions. 

A shift-coupling of random measures (Pi, $i) and (P 2 ,4> 2 ) is a random vector Y 
(possibly on an extension of LI) such that (Pi,dy(4>i)) has the same distribution as 
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(P 2 , $ 2 ). Note that this provides a coupling of the distributions of the two random 
measures such that they are a translated version of each other in every realization. The 
necessary and sufficient condition for the existence of a shift-coupling is proved in 111311 in 
a more general setting. However, the proof in |f!3l is non-constructive (See Theorem |4.9| l. 
Shift-couplings are of special interest in the case that (Pi, $i) is a stationary random 
measure and (P 2 , $ 2 ) is its Palm version. They can be used to construct the Palm version 
of a random measure or to reconstruct the random measure from its Palm version by a 
random translation (see Theorem |4.9[ ). In particular, when $1 is a simple point process, 
Y is a called an extra head scheme for $1 in (5). The name comes from the fact that if 
$1 is the point process on h d formed by the heads in i.i.d. coin tosses (resp. the Poisson 
point process), then by shifting Y to the origin we get the same thing (in distribution) as 
placing a head (resp. point) on the origin. 


3 Motivation: Extra head scheme for a point process 

To get an idea of the definitions in this paper, it is helpful to remind the novel 
algorithm of [3| (Algorithm |3 .1 [ ) which we will generalize. This algorithm is inspired by 
the stable marriage algorithm of Gale and Shapley in bipartite graphs (@) and appears 
to be the first generalization to a continuum setting. The main goal in the paper is to 
construct a translation-invariant balancing allocation between the Lebesgue measure 
and a realization of an ergodic stationary simple point process in R d . 

Let S be a discrete subset of R rf . By sites and centers we mean the points of R d and 
H respectively. The idea of the following algorithm is that each site and each center 
prefer to be allocated as close as possible. Forget the sites that are equidistant from two 
or more centers. 


Algorithm 3.1. For each natural number n, stage n consists of the following two parts: 


(a) Each site x applies to the closest center to x which has not rejected x at any earlier 

stage. 

(b) For each center £, Let A (which depends on n and £) be the set of sites which 

applied to £ in the previous step. Let B be the smallest ball centered at £ such that 
C{A nB)>l. Then £ rejects the sites in A\B. 


Easily seen , each site x either is rejected by all centers or for some center x applies to 
£ and is never rejected for sufficiently large n. In the first case let t(x) = 00 (where 00 is 
treated as a single point added to the space whose distance to every other point is 00 ) 
and in the second case let r(x) := £. 


It is proved in ||3j that the function r defined in this algorithm is stable in a sense that 
is a generalization of stable matchings, which is, roughly speaking, there are no sites 
x\ and x 2 such that \x\ — r(x 2 )| < min{|xi — r(xi)|, |x 2 — r(x 2 )|} (see Definition 4.32 for 
the exact definition). An interesting theorem in |0 is that if T is an ergodic simple 
point process in R d with intensity 1, then almost surely the (random) function t defined 
in Algorithm 3.1 for the discrete set defines a balancing allocation between the 
Lebesgue measure and 4); i.e. almost surely r is defined for almost every site and 
£(t - 1 (£)) = 1 for all centers £. 

One can try to apply Algorithm 1 3. 1| to find an allocation for two arbitrary measures y 
and ip (instead of the Lebesgue measure and a counting measure). For example, let them 
be the Lebesgue measure restricted on (- 00 ,0] and [0, 00 ) respectively. The sites and the 
centers are the points in the supports of ip and ip respectively. In the first step, all sites 
apply to center 0, but the next steps are vague. Same issue may happen in the stationary 
case. If we consider the allocation as a transport kernel T{x, j := S T ( X ), our idea is to 
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force the infinitesimal mass at every point to be spreaded instead of being transported 
to a single point. To do so. We will force an upper bound on the measure T(x, •). Since 
x wants to apply to the support of ip, a natural upper bound is T{x,B) < ip{B). This 
justifies the choice of the name 'constrained' (see definitions |4 .1 1 and |4 .2 [ I. The above 
example will be addressed in Example |6.8[ 


4 Definitions and Main Results 


4.1 Constrained Transport Kernels and Constrained densities 

Let ip and ip be given locally finite non-negative Borel measures on R rf . 

Definition 4.1. A weighted transport kernel T is called ip-constrained ifT(x, •) <ip for 
all x £ R d . 


As we will see in Subsection 4.6 a (ip, ^-balancing ^-constrained transport kernel 
cannot be an allocation except maybe when ip is discrete. 

Since T(x, •) is absolutely continuous w.r.t. ip for any x, we will work with its Radon- 
Nykodim derivative. The following definition and Remark |4.3| establish the setup in a 
measurable way. 


Definition 4.2. A non-negative measurable function f(x,£) on R d x R d is called a sub 
balancing density (given ip and ip) if 


[ /OgOVK^) < l, 

Vx £ R d 

JB. d 


[ f(x,£)ip(dx) < 1, 

)R d 

V£ £ R d . 


We call f balancing if equality happens for ip-a.e. x and ip-a.e. £. Furthermore, f is 
called constrained if it is sub-balancing and f(x, £) < 1 for every (x, £) £ R d x R d . 

Remark 4.3. A constrained density / defines a ^-constrained weighted transport kernel 
via 

T(x,B) := [ f{x,t)iP{dt). (4.1) 

JB 

/(x,£) can be thought as the infinitesimal mass going from x to £. If / is balancing, 
then T is a (ip, ^-balancing transport kernel. Conversely, If T is a (ip, ^-balancing ip- 
constrained transport kernel, then (a suitable version of) the Radon-Nykodim derivative 
of the measure on R d x R d defined as (A x B) i-> J a T(x, B)ip(dx), with respect to ip ® ip, 
is a balancing constrained density. Hence, we might think of / as the density of T w.r.t. 
ip. Similar correspondence holds in the random case. 

Note that the definition of constrained densities depends on both ip and ip, but we 
don't use a prefix for simplicity. The notion of balancing constrained densities is a special 
case of capacity constrained transport kernels defined in OJ. ft can be interpreted as 
there is a transportation capacity between any site and any center. More explicitly, the 
total mass going from a set A of sites to a set B of centers is at most ip(A)ip(B). 


4.2 Construction of a Shift-Coupling 

Let ip, ip £ M be given. Since we have different procedures regarding ip and ip, it 
is helpful to use the names sites and centers for points in the supports of tp and ip 
respectively, following the terminology of Imagine we have two copies of R d , in one 
copy, we have sites and measure ip and in the other copy we have centers and measure 
ip. Nevertheless, we will measure the distance between a site and a center as if they are 
on the same space. We use roman letters for naming sites and Greek letters for naming 
centers. 


Page 6 J25| 




Stable transports between stationary random measures 


Here is an overview of Algorithm |4.4| The algorithm consists of infinitely many stages 
and each stage has two steps. At stage n, each site x 0 applies to the closest possible 
centers with weight A n (x 0 , •) : R d —> [0,1] (tries to construct f{x 0 , •)). Then each center 
Co rejects some of the weights applied to Co if it has reached its capacity. The amount of 
rejection is denoted by C 0 )- Note that even if R„jx 0 , C 0 ) > 0, x 0 will still apply to £0 

at all later stages. The functions A n and R n will be non-decreasing with respect to n. 
A heuristic for choosing A n and R n in the algorithm is that the sites prefer to apply to 
centers which are as close as possible, in a greedy manner. Similarly, the centers prefer 
to reject no portion of the applications of the sites which are as close as possible. 

For an illustration of Algorithm |4.4| see Example |6.2[ Example |6. 1 1 shows that this 
algorithm generalizes the algorithm in (3j. The name site-optimal is justified in Corol¬ 
lary [4?22] 

Algorithm 4.4 (site-optimal Gale-Shapley algorithm). Given measures and d on R d , 
let the rejection function be zero at the beginning; i.e. R 0 (x,^) := 0 for all ( x , C) G R d x R d . 
For each natural number n, stage n consists of the following steps: 

(i) For each site x 0 , define its application radius at stage n as 


i(®o) := sup la: (1 - £)) ip{d£) < 1 > . 

[ JB a (x 0) ) 


Define the ?z-th application function as 


(4.2) 


A n {x o, C) 


1 |®o - Cl < a n (x 0 ), 

cR n - i(x 0 ,£) + (1 - c) |x 0 — £| = a n (x 0 ), 
0 l^o - Cl > a n (x 0 ), 


where c = c n (xo) is the constant in [0,1] such that we have 

/ (A n (x 0 ,O - Rn-i(xo,0) fp(d£) = 1 ifa n (x 0 )<oo 


(4.3) 


Jn d 

and we let c = 1 (or any arbitrary constant) if a n (x 0 ) = 00 or t/>(< 9R Qji ( Xo ) (x 0 )) = 0. 
We say x 0 applies to the centers with weight A n (x 0 , •) at stage n. 

(ii) For each center Co, define its rejection radius at stage n as 


r n (€o) := sup 


' B r (% 0 ) 


A n (x,£ 0 )ip(dx) < 1 


Define the n-th rejection function as 


Rn £ 0 ) 


0 ]a: - Co| < r„(Co), 

c'A n (x,£ 0 ) \x - Co| = r„(Co), 
A n (x,£, 0 ) \x - Co| > r n (Co), 


where d = d n (£ 0 ) is the constant in [0,1] such that we have 


(4.4) 



(A n (x, Co) - Rn(x, Co)) <p{dx) 


1 ifr n (£o) < OO 


and we let d = 0 if r n ( Co) =00 or ^(dB rn ^ 0 ) (Co)) = 0. We say Co rejects the 
application weights according to R n (-, Co) at stage n. 


It is shown in the proof of Lemma 4.5 that c n (xo) and d n ((jo), which are defined in the 
algorithm, exist and are well defined. 
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Lemma 4.5. In Algorithm \4.4\ A„ , R n and a n are non-decreasing with respect to n, r n is 
non-increasing. Moreover, they depend on p and uj in a measurable and How-adapted 
manner. 

This lemma allows us to provide the following definitions. 

Definition 4.6. In Algorithm \4.4\ the (final) application radius, rejection radius, applica¬ 
tion function, rejection function and the site-optimal density are defined as follows 


a(x) 

:= lim a„(x), 

n—>oo 


:= lim r ■„(£), 

n—>o o 

A(x,{) 

:= lim A„(:r,£), 

n—>■ oo 

R(x,0 

:= lim R n (x,^), 

n—t oo 

fs(x ,0 

:= A(x,£) - R{x,0 


Definition 4.7. The center-optimal Gale-Shapley algorithm for p and ip is just Algo¬ 
rithm^^ for ip and p; i.e. we swap the roles of sites and centers. The center-optimal 
density f c is defined similar to Definition |4.6| using the center-optimal Gale-Shapley 
algorithm. 


Theorem 4.8. Let $ and 4> be ergodic random measures on with positive and finite 
intensities. If the intensities are equal, then the site-optimal density which is constructed 
by Algorithm \ 4. 4] is almost surely a balancing constrained density. Therefore, it almost 
surely gives a flow-adapted (4?^, 4^,) -balancing transport kernel via ( |4. 1 [ > . 

According to Theorem |4.8| Algorithm |4.4| answers a question in |[8l asking for a 
constructive algorithm to find a flow-adapted balancing transport kernel, although the 
existence of such a transport kernel is proved to be equivalent to the equality of the 
sample intensities in (8j, which is granted here by ergodicity. 

As a result of Theorem |4.8| one can construct a shift-coupling between an ergodic 
random measure and its Palm version as in the following theorem. The key tool for this 
construction is the following well-known theorem in the literature. 


Theorem 4.9 (Shift-coupling). Let 4/ be an ergodic random measure on R d with positive 
and finite intensity A>i. and let 4> = X^Cd- Let F = F u (x, £) be a flow-adapted function 
which is almost surely a balancing density for 4> w and 4> w ; e.g. the site-optimal density 
for 4> w and 4c,. 


(i) IfY is a random vector such that its conditional distribution given 4/ is F((). •)'!', 
then Y gives a shift-coupling of 4/ and its Palm version; i.e. Oy T has the same 
distribution as the Palm version of T. 


(ii) On the probability space (fi, T, ll\[,), ifY is a random vector such that its conditional 
distribution given 4/ is F(-, 0)$, then the random measure (P*, Oy'to) has the same 
distribution as (P, 4/). In words, (fy T is a reconstruction of 4/ from its Palm version. 

Theorem |4.9| is a direct implication of Corollary 4.7 in O (see also Theorem 16 in [0 
for point processes). Nevertheless, we state it for the purpose of this paper. In part Q, 
we call the random vector Y an extra head scheme for 4/ with the terminology of @. 
One can also replace F(0, •)4' by T{ 0, •), where T is a flow-adapted transport kernel 
which is almost surely (4>, 4/)-balancing. If moreover T is a balancing allocation that is a 
function of 4', provided that it exists, then Y will be a deterministic vector conditional 
on 4>. This is called a non-randomized extra head scheme (note that to construct 
Y in the general case, one may need extra randomness; i.e. to extend the probability 
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space) and its existence is proved in the case when is an ergodic simple point process 
in |3ll by providing a construction. Algorithm 1 4.4 1 is a generalization of that construction. 
Although it does not provide a non-randomized extra head scheme in Theorem |4.9| it has 
the property that the distribution of Y conditional on 'h depends only on the realization 
of 'll. We don't know whether non-randomized extra head schemes always exist for 
general random measures or not, as stated in Open Problem |4.39[ However, we show 
that a flow-adapted balancing allocation may not exist for two general random measures, 
as shown in examples |6.6| and |6.7| Some results and open problems about this problem 
are provided in Subsection |4.6| 

The claim of Theorem |4.9| is similar to the inversion formula ( 111111 . Satz 2.4. See 
also (8]|, (2.6)). This formula recovers the distribution of T from its Palm distribution. 
Here, the difference is that the balancing property of F ensures that the recovery can 
be done using a shift-coupling as defined in the theorem. 


Remark 4.10. In Algorithm |4.4| one could define A n (x o, •) on dB, 
on dB rri (Co) in other ways. If this is done such that Lemma 
then all of our results remain valid. 


4.5 


„(xo) (zo) an d R n (*> Co) 

and Remark!! 


5.1 


hold. 


Remark 4.11. Being stationary is crucial in Theorem |4.14| As an example. Example |6.3| 
shows cases in which the site-optimal density is not balancing although the (non-random) 
measures in the example may have equal spatial intensities; i.e. y>([—r, r]) ~ ip([—r,r]) 
as r — > oo (see also CD). As mentioned in it seems difficult to give a sufficient 
condition for deterministic discrete sets (and therefore, for measures) to ensure that the 
site-optimal density is balancing. 


4.3 Stability of Constrained Densities 

Definition 4.12. Let ip,ip £ M be given and f be a given constrained density. We say 
that a site * 0 is f -exhausted if 

I f(xo,£,)il>(d£) = 1 

JR d 

and f -unexhausted otherwise. Similarly a center Co is f -sated if 

[ f{x,Zo)<p(dx) = l 

JR d 

and f -unsated otherwise. We say that a site x 0 , f -desires a center Co if f{x o,£o) < 1 
and either x 0 is f-unexhausted or 

3£i G R d : \x 0 - Ci| > |o:o -Co| and /(* 0 ,Ci) > 0. 

Similarly, we say( 0 , f -desires x 0 if f{x o,£o) < 1 and either Co is f-unsated or 

3*i e R d : |*i - Col > \%o - Col and /(*i,Co) > 0. 

We drop the prefix '/-' when there is no confusion. A constrained density f is called 
stable if there is no (* 0 , Co) G R d x R d such that * 0 desires Co and Co desires * 0 . 

In this definition, each sites prefers the centers according to Euclidean distance 
and vice-versa. See examples |6.3| and |6.4| in Section[6]for examples of a stable and an 
unstable constrained density. 

Theorem 4.13 (Stability). The site-optimal density is a stable constrained density and 
depends on ip and uj in a measurable and flow-adapted manner. 

Theorem 4.14. Let $ and T be ergodic random measures on R d with positive and Unite 
intensities A$ and Ay. Let F = F ul (x, C) be a flow-adapted function which is almost surely 
a stable constrained density for <1>^. and Then almost surely 
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(i) if A$ = then F u is (<& u , 4> u )-balancing; i.e. the set of unexhausted sites has zero 

-measure and the set of unsated centers has zero -measure. 

(ii) if A,[, < A®, then there is no unexhausted site but the set of unsated centers has an 
infinite -measure. 

(iii) if A,j, > A®, then there is no unsated center but the set of unexhausted sites has an 
infinite -measure. 


Proposition |4 .2 0| quantifies how far the centers (resp. sites) are from being sated 
(resp. exhausted) in average in the second (resp. third) case of Theorem |4.14| 

Remark 4.15. Theorems 4.8 4.9 and 4.14 can be extended to stationary non-ergodic 
cases. The necessary and sufficient condition on <f> and T is the equality of the sample 
intensities linp—xx) 4>(B r (0 ))/C(B r (0)) and lim^oo 4t(B r (0 ))/C(B r (0)). This condition 
is proved to be necessary and sufficient for the existence of a flow-adapted balancing 
transport kernel in Theorem 5.1 in 10. To do so, we can slightly change the proofs and 
use conditional expectation with respect to the invariant sigma-filed (using the ideas 
in (8]|) or obtain the general versions simply by applying the ergodic decomposition 
theorem. 


4.4 Other Properties of Stable Constrained densities 

In this subsection we study monotonicity, optimality, uniqueness and boundedness 
of territories of stable constrained densities. In general, uniqueness is not granted in 
the deterministic case (see H3j for a counter example). The manner of choosing between 
equidistant points in Algorithm |4.4| (see Remark [4.10 1 is another obstacle for uniqueness; 
e.g. when <p and ip are measures on Z d . Assumption 4.16 gives a sufficient condition for 


uniqueness of the choice in Remark [440] for almost all points. However, we will prove 
uniqueness of stable constrained densities only in the stationary case in Theorem |4.24| 
It can be seen that the first condition in Assumption |4.16] means that, for ip-a.e. site x, 
the boundary of no ball centered at x can be partitioned in two disjoint sets with positive 


ip-measure. This assumption is not difficult to satisfy as shown by Proposition 4.17 
Assumption 4.16. For ip-a.e. site x and ip-a.e. center £, we have 


Vr >0, 3s € dB r (x) : ip(dB r (x) \{s}) = 0, (4.5) 

Vr >0, 3s G dB r (£) : tp(dB r (£) \{s}) = 0. (4.6) 


Proposition 4.17. If at least one of ip and ip assigns zero to all spheres and all affine 
hyperplanes of R d ; e.g. is absolutely continuous with respect to the Lebesgue measure, 
then Assumption\4. 16\ holds. 


The following auxiliary functions measure how far the mass is transported from sites 
and to centers, given a sub-balancing density /. 


Definition 4.18. Let f be a sub-balancing density as in Definition 4.2 Fort e [0, oo] 
define 


x(f,t ) := [ f(x,£)l\ x -e\< t ip(d(;), 

JR d 

■■= I f (x, £)l\ x -£\< t <p(dx). 

JR d 

Remark 4.19. If /(ar,£) < 1 and ip x (f, \x - £|) < 1, then x desires £. Similarly, if f(x ,£) < 
1 and \x — £|) < 1, then £ desires x. Note that the converse is not true; i.e. x may 
desire £ even if ip x (f, \x - £|) = 1 since f(x, •) can be positive on a i/)-null set that contains 
a center farther than £ to x. 
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Proposition 4.20. In the setting of Theorem 4.14 one has 


lim 

r—> oo 


lim 

r—too 


1 

W (0)) 

i 

m 



T' x (F,oo)$(dx) 


'BA 0 ) 


&(F, oo)tl'(d^) 


r | . 

mm{l, -— }, a.s. 

Acj> 

r | , 

mm{l, -—}, a.s. 

Avp 


Note that the left hand sides in the above equations are the spatial averages of 
dl x (F, oo) and &(F, oo). This is a measure of how far the sites and the centers are from 
being satisfied. 


Theorem 4.21 (Monotonicity). Let (p,ip) and (p, v) he two pairs of measures such that 
A > T>, v < ip and the pair (p, ip) satisfies Assumption 4.16 Let f be an arbitrary stable 
constrained density for (p, v) and consider the site-optimal density f s for (ip, ip) together 
with the functions in Definition \4.6l 


(i) We have 


f + R< 1, (pgip)-a.e. 


(ii) For (p g ip)-a.e. (x,pP), if x fully applies to for example if \x — < a(x), then 

fM < f s (x,0- 

(Hi) For p-a.e. site x we have 


Vf€[0,oo]. 

(iv) Forip-a.e. center £ we have 

VfG[0,oo]. 


Intuitively, it means that when there are less sites and more centers, the situation is 
better for sites and worse for centers. The following corollary is immediately obtained 
from Theorem |4.21[ Note that the inequalities that contain f c are equivalent to the ones 
that contain f s as seen by swapping the roles of the sites and the centers. 


Corollary 4.22 (Optimality). Suppose p and ip satisfy Assumption 4.16 and let f be an 
arbitrary stable constrained density for ip and ip. Let f s and f c be the site-optimal and 
the center-optimal densities for the same measures. 


(i) For p-a.e. site x we have 


lpx(fs,t) > 1p x (f,t ) > 1p x (fc,t), Vf G [0, oo]. 


(ii) Forip-a.e. center £ we have 

<P*(fc,t) > pHf,t) > p £ (f s ,t), Vf G [0, oo]. 

In words, among all stable constrained densities, the site-optimal density is the 
best for sites and the worst for centers. This justifies the names site-optimal and 
center-optimal for the Gale-Shapley algorithm. 

As mentioned in @1, it seems difficult to express a simple condition in terms of p 
and ip that ensures uniqueness of stable constrained densities. But As sumption |4 .1 6] 
is enough in the stationary case, as shown in Theorem |4.24| The key for proving the 
theorem is the following proposition. 
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Proposition 4.23. With the assumptions of Corollary\4.22\ if for ip-a.e. site x we have 


ipx{fs,t)=il> x (fc,t), Vt£[0,oo], 


then there is a (ip® ip)-a.e. unique stable constrained density for p and ip. 

Here, by (ip <S> ip)-a.e. unique, we mean that any two stable constrained densities are 
identical except on a set of zero (ip <g> ^-measure. 

Theorem 4.24 (Uniqueness). Let$ and be stationary random measures on R d with 
positive finite intensities that satisfy Assumption \4. 1 6 | almost surely. Almost surely, any 
two stable constrained densities for and '10 agree on (4> w 0 'l 0 ,)-a.a. points. 

Here, we mean that there is an event with probability one such that uniqueness holds 
on that event. 

Given a stable constrained density /, the territory of a site x is the set {£ £ R d : 
/(x,£) > 0}. Similarly, the territory of a center £ is the set {a; £ R d : /(x,£) > 0}. 

Theorem 4.25 (Boundedness). Let$ and be stationary random measures on R d which 
are almost surely non-zero and have finite intensities. LetF = F w (x,£) be a flow-adapted 
stable constrained density for (4>, T). Almost surely we have 

(i) 40-a.a. sites and 4> w -a.a. centers have bounded territories. 

(ii) The union of the territories of the sites (resp. centers) in a bounded set, has finite 
Hi-measure (resp. finite 4>- measure). 

4.5 Voronoi Transport Kernel Corresponding to a Measure 

In this subsection, we generalize the notion of Voronoi tessellation and define it for a 
measure. We will use it only for proving Theorem |4. 2 5| here, but it can be interesting in 
its own. 

Definition 4.26. Let ip be a measure on R d such thatip( R d ) > 1. For xg £ R d define 


s(x q) := sup{s £ R : ip(B s (x 0 )) < !}• 


Forx 0 ,£ £ R d define 



v(x 0 ,€) := < c, ]aj 0 - Cl = s(x 0 ), 
[ 0, |x 0 - Cl > s(x 0 ), 


where c is the constant in [0,1] such that J v(xq, Oi’idO = 1. We let c = 1 if s(x 0 ) = oo 
or ip(dB s ( Xo ) (x ())) = 0. The Voronoi density and the Voronoi transport kernel with 
respect to ip are the function v and the transport kernel V(x, •) := v(x, ■)ip respectively. 
The Voronoi territory of center f with respect to ip is the set {x £ R. f/ : v(x, p) > 0}. 

Note that the Voronoi density is the same as the function A\ in Algorithm |4.4| 
Remark 4.27. The condition ip(R d ) > 1 ensures that the Voronoi transport kernel is a 
non-weighted transport kernel. If ip is the counting measure on a discrete set 2, then the 
Voronoi territories with respect to ip give the usual Voronoi tessellation of 2. Therefore, 
the notion of Voronoi transport kernel generalizes the notion of Voronoi tessellation. 

Remark 4.28. It is easy to see that x is in the Voronoi territory of £ if and only if 
ip{B? £i (x)) < 1 and if equality happens, then ip(dB\ x _^ (x)) = 0. 

Lemma 4.29. The Voronoi territory of £ is star-shaped with center £ but not necessarily 
convex or closed. Moreover, it is a closed polyhedral region in R d provided that ip has a 
discrete support. 
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Proposition 4.30. If every half-space has infinite V’- measure, then all Voronoi territories 
with respect to ip are bounded. In particular; if 'I' is a stationary random measure which 
is almost surely non-zero, then almost surely all Voronoi territories with respect to are 
bounded. 


4.6 On Allocations 


As explained in the introduction, a problem of special interest is constructing flow- 
adapted balancing allocations. But since a ^-constrained weighted transport kernel T 
satisfies T < ip (as in Definition 4.11, a balancing ^-constrained transport kernel cannot 
be an allocation except maybe when ip is a discrete measure. The reason is, if B is 
a set of centers such that 0 < ip(B) < 1 (which exists in the non-discrete case), then 
on a set of sites with positive measure (under p) we have T(x , B) > 0. On the other 
hand, T(x, B) < 1 and so T(x, •) is not a Dirac measure. However, we will see that under 
some conditions in the case when is a counting measure, the site-optimal density is 
guarantied to give an allocation. 

Motivating from f3j and Definition [47X2] we define 


Definition 4.31. Given p,ip £ M, an allocation is a measurable function r : D —> 
R d U {oo}, where D C R d is a measurable set such that p(D c ) = 0 and oo is treated as 
a single point added to the space. This can be regarded a weighted transport kernel 
that transports the mass at s to a single point r(s); i.e. define T(s, •) to be < 5 r ( s ) whenever 
s ^ DU t~ 1 ( oo) and zero otherwise. This allocation is (ip, ip)-balancing if r ^ oo, ip-a.e. 
andpir- 1 ^)) =ip(-). 

Definition 4.32. Given measures ip and ip on R d , let r be an allocation defined on a 
domain D C supp(ip) with p(D c ) = 0. We say r is sub-balancing when p(r~ 1 (B)) < 
’ip(B) for all B e Q; i.e. t*(p\d) < ip. Given a sub-balancing allocation r, a site xi, € 1) 
is exhausted when r( xo) oo and unexhausted when t(xq) = oo. A center Co is sated 
when it is not in the support of the measure ip - r*(p\ D ). we say a site x 0 desires a 
center £ 0 when either x 0 is unexhausted or \x 0 — t(x 0 )| < \x 0 - Col- We say £ 0 desires 
x 0 if t(x o) 7 ^ Co and either Co is unsated or there is a sites x\ s.th. t(x i) = Co and 
\xi - Co| > |xo — Col- Wfe say that a sub-balancing allocation is stable if there is no pair 
(x 0 ,Co) such that both desire each other. 


A similar notion of one-sided stable allocations in dimension one is also defined in |[9]|. 
The authors construct a balancing stable allocation (in that sense) between two ergodic 
stationary and mutually singular diffuse (i.e. without atom) random measures on R. 


Proposition 4.33. In the following cases the site-optimal density is {0, l}-valued on 
(<p ® ip)-a.e. points. 


(i) when both p and ip are absolutely continuous w.r.t. C, 


(ii) when one of p and ip is absolutely continuous w.r.t. C and the other is a counting 
measure. 


(Hi) when both p and ip are counting measures and Assumption 4.16 holds. 

We skip the proof of Proposition [4B3] For a proof, it is easy to use induction to show 
that the functions in Algorithm [4A] are {0, l}-valued a.e. 

Remark 4.34. In the case that each one of p and ip is either a diffuse (i.e. without 
atom) measure or a counting measure (not necessarily satisfying Assumption |4.16[ l. 
Algorithm 4.4 can be slightly changed to find a {0, l}-valued stable constrained den¬ 
sity. It is enough to re-define the functions A n (x 0 , •) and i?„(-,Co) on dB a t Xo ) (x 0 ) and 
dB rn(io) (Co) for example using the lexicographic order on the boundary of the balls (see 
Remark |4. 10 P . 
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Proposition 4.35. A constrained density which is {0, l}-valued on (p <g> ip)-a.e. points 
gives an allocation provided that ip is a counting measure. Therefore, by Remark 4.34 
for the existence (and construction) of stable allocations it is enough that 

(i) p is diffuse and ip is a counting measure, 

(ii) or both p and ip are counting measures. 

Here, by 'gives an allocation' we mean that the weighted transport kernel given by 
Remark l4. 31 coincides with an allocation on almost all sites. 

Remark 4.36. When ip is a measure with discrete support ip = Xiti and p 
is diffuse, we can slightly change Algorithm |4.4| to obtain a stable allocation. We 
limit ourselves to ^-constrained transport kernels, where ip = YlrLi %, is the counting 
measure with the same support as ip. Equivalently, f{x,£i) < d- for defining constrained 
densities. We then change the definition of A n and a n in Algorithm |4.4| in a similar 
manner. We similarly get the following proposition. 

Proposition 4.37. Suppose *!> is a diffuse random measure and 'k is a random measure 
whose support is a.s. discrete in the setting of Theorem \4.8\ One can construct an 
allocation which is almost surely stable and (<i>. 4/) -balancing. 

We finish this section with two open problems. 

Open Problem 4.38. Does there exist a stable allocation for any two measures p and 
ip on R d provided that p is absolutely continuous with respect to the Lebesgue measure? 

If a stable allocation exists and can be chosen as a measurable function of (p, ip), then 
it leads to a solution of the following question in a way similar to theorems |4.8| and |4.9| 

Open Problem 4.39. Does there exists a non-randomized extra head scheme for any 
ergodic random measure on R d with a positive and finite intensity? 

5 Proofs 

Proof of Lemma IZ 71 Assume P[4>(V) < oo] = 26 > 0 an let e > 0 be arbitrary. This 
implies that there exists R > 0 such that P[4 >(V"\.Br( 0)) < e] > 5. For a given r > 0, 
let B r C V\Br( 0) be a ball with radius r, which exists by the assumption on V. By 
stationarity we get 


P[$(B r (0)) < e] = P[$(B r ) < e] > P[$(V\B fl (0)) < e] > <5. 


By letting r —> oo we get P[4>(R d ) < e] > <5. Since e is arbitrary, this gives P[$(R d ) = 0] > 
5, a contradiction. 


O 


Remark 5.1. By the definition of A n we have 



and if a n ( xq) < oo, equality holds. Moreover, if ip has infinite total mass, then one can 
prove by induction that a n (-) < oo for all n. Similarly, by the definition of R n we have 



and if r n (£ 0 ) < oo, equality holds. Moreover, if equality holds at some stage n, then it 
holds at all later stages; i.e. £ 0 is sated at stage n afterwards. 
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Proof of Lemma [431 It is clear that R± > Rq. For n > 2, assume R n -\ > R n ~ 2 - We will 
conclude a n > a n ~\, A n > A n -i, r n < r n _i and which proves the claim by 

induction. By the same induction, it is easily proved that the functions depend on ip and 
ip in a measurable and flow-adapted manner. 

Since 1 — R n -i < 1 — R n - 2 , H4.2| l gives a n (x 0 ) > a n -i{x 0 ). If we have a n (x 0 ) > a n -i(x 0 ) 
or a n (x 0 ) = oo, then it is clear from the definition of A n that A n (x 0 , •) > A n _i(ar 0) •)• Now, 
suppose a n (x 0 ) = a n _i(x 0 ) =:t< oo. By ( |4.2p and H4.3| l we get 

Cn[X0) f gB (l-R n -i(x 0 ,O)m) ’ 

where B = B t (x 0 ) and jj = 0 by convention. This equation and the fact that R n -i > Rn -2 
implies c n {x 0 ) < c„_i(a;o) and hence A n (x 0 ,-) > A n - Given A n > A n -i, H4.4} 
implies r„(( 0 ) < r n -i(Co)- The proof of the fact that r n < r n - 1 implies R n > R n _ 1 is 
completely similar to the above proof. 

□ 


Lemma 5.2. In the site-optimal density, if a site xo desires a center Co, thenA n (x 0 ,Co) = 1 
for sufficiently large n. Similarly; if£ 0 desires x 0 then R n (xo,£o) = 0 for all n. 

Proof. Suppose x 0 desire £ 0 - If there is a center £i such that |^i — x 0 | > |^o — x 0 \ and 
f(x o,Ci) > 0- then x 0 has applied to ^ and so it has applied to all points closer than 
with weight 1. Hence xq has applied to £o with weight 1. In the other case, xq is 
unexhausted. It is enough to prove that a n (x 0 ) > |xo - ^o| for sufficiently large n. If this 
is not true, a n (x 0 ) < oo for all n and Remark [5H] gives 



{A n (: r 0 ,£) - Rn-i(xo,€))ip(d€) = 1. 


Since a n ( x 0 ) < |®o - Co|, A n (x 0 ,-) - is bounded by ls |a , 0 _ 4ol (* 0 ), whic h is 

integrable with respect to ip due to locally finiteness of ip. So, Lebesgue's dominated 
convergence theorem gives f Rd f s (x 0 ,0' l l J (d() = 1; i.e. x 0 is exhausted, a contradiction. 

Now, suppose that £ 0 desires x 0 . If £ 0 is unsated, then it has not rejected any weights. 
So suppose there is a site x\ such that \x\ — £ 0 | > \ x o — ^ol an< i f{% i,^o) > 0- Therefore 
7Vi(£o) > \x\ - Col for all n, or else Co would fully reject x\ at all stages after stage n since 
the rejection radius is non-increasing. Thus r„(Co) > |a; 0 — Col for all n and so Co has not 
rejected any application weight of x 0 . □ 


Proof ofTheorem \4.13l Let f s be the site-optimal density. By Lemma 4.5 we get that f s 
depends on c p and to in a measurable and flow-adapted manner. We have 


f s = lim A n - R n - 1 = lim A n - R n . 

n—> oo n—>■ oo 


Moreover, A n > R n -\ and A n > R n ■ By Remark 5.1 and Fato's lemma, it follows that 


f s is a sub-balancing density. Also, f s is constrained since f s < 1. Now, suppose that 
f s is unstable. So we can find a site x 0 and a center Co that desire each other. By 


Lemma 5.2 we get that A„(a:o,Co) = 1 and R n (x o,Co) = 0 for sufficiently large n. So 


fs( x o,£o) = 1 - 0 = 1, a contradiction. 


□ 


Lemma 5.3. Let f be any stable constrained density for p and ip as in Definition \4.12\ 
Let X' be the set of unexhausted sites and S' be the set of unsated centers. If X' ^ 0, 
then ip(B,') < 1 and if S' ^ 0, then p(X') < 1. In particular, we have either <p(X') < 1 or 
V-(s') < l. 
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Proof. First, suppose X' ^ 0 and ip(E') > 1. Let a- 0 € X'. Since x$ is unexhausted and 
V'(S') > 1, we find a point £ 0 £ S' such that f(x 0 ,£, o) < 1- Now (x 0 ,fy) is an unstable 
pair since x 0 is unexhausted and is unsated, a contradiction. Similarly, if E' / 0 we 
conclude that p(X 0 ) < 1, which completes the proof. 

□ 


Note that it is possible that both X' and 3' in Lemma [5,3| have positive measure, as 
shown in Example |6.2[ 


Lemma 5.4. Let $ and if be stationary random measures on lR d satisfying the assump¬ 
tions ofTheorem\4.14\ For any flow-adapted sub-balancing density F we have 


A$E$ [*o (F,t)} = A*E* [<h°(F,t)] . 


(5.1) 


Intuitively, this means that the average mass that is transported to a typical center 
from the sites of distance at most t is equal to the mass that is transported from a typical 
site to the centers of distance at most t. This is a version of mass transport principles. A 
more general equation can be found in If8l . 


Proof. Let 7T(w, s) := F UI (Q, s)l| s | <t . By Neveu's exchange formula ( |2.4| > we get 


A$E$ [T'o(E,f)] 



[ Fg o (0,s)l\ s \ <t 'L{ds) 
,jR d 


A^E^ 

/ Fg s (0, —s)l|_ s | <t 4 , (ds) 
J R d 

Av/aE^ 

/ Fe o (s,0)l\ s \<t^(ds) 

JR d 



A*E* [$°(F,f)] . 


□ 

Proof of Theorem 14. 1 41 Let and ifi be the restrictions of d> and if to the set of 
unexhausted sites and unsated centers respectively. Since F is flow-adapted, $1 and if 1 
are also ergodic random measures. Using Definition |4.18| define 

U X {F) := l-f,(F,oc), 

U^(F) := 1-$£(F, 00 ). 

By Lemma [5~4] for t = 00 we get 

A$E$ [1 - J7 0 (F)] = A^Evt [1 - U°(F)] . 


Thus 

A*E* [J7 0 (F)] - A*E* [17° (F)] = A* - A*. 
By \2.2) , if B is a cube in R d , then we have 

1 


Cd(B) 

1 

C d (B) 

1 

C d (B) 


IJB 


U 0 (F ea )$(ds) 

U s (F)<h{ds) 

U s (F)^ 1 (ds) 


(5.2) 
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where the last equality is due to the fact that $1 = 1{e/ s (f)>o}3l Since U S (F) is bounded 
and $1 is ergodic, we get that 


E$[(7 0 (F)]=0 o 
E* [C/bOF)] > 0 «• 


E [d>i(f3)] = 0 <£> $i(R d ) = 0, a.s. 
E [<E>i(i?)] > 0 $i(R d ) = oo, a.s. 


Similarly 


E*[£/°(F)]=0 o 
E* [C/°(F)] > 0 


E [^i(R)] = 0 ^ 4-i(R d ) = 0, a.s. 
E [^i(R)] > 0 ^i(R d ) = oo, a.s. 


(5.3) 


(5.4) 


Since F is stable almost surely Lemma 5.3 gives that, almost surely either $i(R d ) < oo 
or \&i(R d ) < oo. By ergodicity either 4>i(R d ) < oo a.s. or 4'i(R d ) < oo a.s. Hence, 
according to ( |5.3[ > and ( |5.4[ >, either E$ [Z7 0 (.F)] = 0 or E^ [(7°(F)] = 0. If we substitute 
this in ( f5~2l >, we get 

' ®*[Uo(F)] = max{0,l-ff}, 

E*[J7°(F)] = max{0,l-^}. 

These equalities, ( |5.3| l and ( |5.4| l complete the proof. □ 


Proof ofTheorem \4.8l By Lemma [475] the site-optimal density for and T gives a flow- 
adapted transport kernel. Moreover, it is a stable density by Theorem |4.13| Now the 
theorem is a direct consequence of part 10 of Theorem |4.14| □ 


Proof ofTheorem \4.9l Let H be any measurable function on M. Since P$ = P, it is a 
direct consequence of Corollary 4.7 in j8j that E^ [if (\D)] = E [If (0y\P)] in case 0 and 
E^ [if(0y\l/)] = E[if(4')] in case 0). This finishes the proof. In fact, case 0 is just 
Example 4.8 in l|8| . 

□ 


Proof of Proposition ^ . 1 7| Suppose ip assigns zero to all spheres and all affine hyper¬ 
planes. So p(dB) = 0 for every ball B. Thus, ( [4761 ) holds for all £. For 0 < k < d, 
by a fc-dimensional sphere, we mean the intersection of the boundary of a ball with a 
non-tangent affine subspace of dimension k + 1. We call a fc-dimensional sphere S bad if 
-0( < S'\{s}) > 0 for every point s. A bad sphere is called minimal if it contains no other bad 
spheres of lower dimension as a subset. If ( |4.5| ) fails for a site x, then x is equidistant 
from the points of a minimal bad sphere. Since the set of points that are equidistant 
from all points of a sphere is a proper affine subspace of R d , it suffices to show that 
there exist only a countable number of minimal bad spheres. 

Let T be the set of atoms of ip, which is countable due to locally finiteness of ip. 
zero-dimensional bad spheres are just pairs of atoms and so they are countable. Also, a 
positive-dimensional minimal bad sphere S contains at most one atom and so ip(S\T) > 0. 
Suppose there are an uncountable number of minimal bad spheres. So, there exist R > 0 
and e > 0 such that there are infinitely many minimal bad spheres §j,i £ IN such that 
Si C B r (0) and ip(Si\T) > e. For i j, since Si D Sj is not a bad sphere, we have 
ip (( Si\T ) n ( Sj\T )) = 0. Therefore ip(B R (0)) > 07 ip(Si\T ) = oo, a contradiction. □ 

Proof of Proposition \4.20\ The claim is a direct consequence of H5.5P and Birkhoff's 
theorem. □ 

Proof of Theorem \4.21\ 0 Suppose the statement is false and let n be the first stage that 
/(x,£) + Rn(x, £) > 1 for a positive (p <g> ■i/ij-measure of pairs (a;,£). By Fubini’s theorem 
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and Assumption 4.16 for p and if) we can find a set 
center £ 0 £ Si, the set 


E/ with V'(-i) > 0 such that for each 


T '■= T to - = { x ' f( x > Co) + Rn ( x , Co) > 1} 


has positive p-measure and 

Vr > 0,3s £ dB r (£ 0 ) : p( dB r (£ 0 ) \{s}) = 0. (5.6) 

By the definition of T, we have R n (x, Co) > 1 — f( x , Co) > 0 for x £ T. So £ 0 has rejected 
some weight from all sites in T at stage n and so £ 0 is sated at that stage. Thus, if we let 

B :=B io :=B rn{ia) (£ 0 ), we have 


/ A n (x, £o) - Rn(x, Co )<p{dx) = 1 
J B 


(5.7) 


and moreover, T is disjoint from the interior of B. Since f Rd f(x, £o)p(da:) < 1 and 
p > <p, < [5771 ) gives 

f [A n -Rn- f](x,€o)<p(dx) > 0. 

J B 


(5.8) 


Lemma 5.5. There is a subset Be C B with positive y- measure such that for all x £ B' £ 
we have 


So 


(a) [A n - R n - f](x,£ 0 ) > 0, 

(b) x is closer than some point in T to £ 0 . 

To prove the lemma, we consider three cases. Note that if p is absolutely continuous 
w.r.t. the Lebesgue measure, then only the first case happens. 

Case 1. Suppose p(T(ldB) = 0 and thus p(T\B) > 0. Since p(T) > 0 and /(x, Co) > 
1 - Rn{x, £o) > 0 for all x £ T, we have f T f(x, £o)p,(dx) > 0. Thus, by the assumption of 
this case, 

[ f{ x ,to)v(dx) < 1 - f f(x,€ Q )p,(dx) < 1. 

J B JT 

Therefore, the inequality in ( ]5.8[ ) is strict. Thus, the integrand is positive on a set with 
non-zero p-measure, which is the desired set. 

Case 2. Suppose p(T n dB) > 0 but <p(T n dB ) = 0. By ( |5.6P we get p(dB\T) = 0 and 
thus ip(dB\T) = 0. So tp(dB) = 0 and we can replace B by B° in ( |5.7D and ( |5.8| l. The 
rest of the argument is similar to the previous case since B° n T = 0. 

Case 3. Suppose <p(T D dB) > 0. Since / + R n > 1 on T x {£o}> we have 


[ [A n - R n - f](x, Co )<p{dx) < 0. 

J TndB 


On the other hand, ( f5~6] > gives ip(dB\T) = 0. Now (578]) gives 


[ [A n - R n - f](x,£o)<p(dx) > 0. 
J B° 


So the integrand is positive on a set with positive y-measure, which is the desired set. 
This completes the proof of Lemma [575] 
part © of Lemma [ 


For x £ B'c o , 


5.5 


implies f(x, £ 0 ) < 1- By the definition of T we get 


/(■, Co) > 0 on T. So £ 0 , /-desires x because of part (bj of Lemma 5.5 i.e. it desires all 


sites in B'^ q . 
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By Fubini's theorem for the set {(x,£) : £ £ Ei,x £ B £} and Assumption 4.16 for p 
and if we get that there is a site x 0 such that the set 


{£ £ si : Xq £ B'^j 


has positive /’-measure and 

Vr > 0,3s £ dB r (x 0 ) : if(dB r (x 0 ) \{s}) = 0. (5.9) 

Note that our construction of B^ o is given in terms of some inequalities. So the above set 
is measurable and Fubini's theorem is valid. Let Bi be the smallest closed ball centered 
at xq (with possibly infinite radius) that contains E 2 . Part © of Lemma [53] implies 

1 - f?„_i(a" 0 ,£) > [A n -R n - i]Oo,£) (5.10) 

> [An - f?„](x 0 ,£) 

> fix o,£)>0, V£ £ S 2 . 

As a result, x 0 has applied to all centers in S 2 at stage n. Therefore, A n {x 0 , •) = 1 on 
BI. Furthermore, by the choice of n we could choose x 0 such that if-a.e. we have 

f(x o, •) + R n _i(x 0 , •) < 1. Thus 


— _n„_ipxo) ■) ' 
If ip(E 2 (~l Bf) >0, ( |5. 1 0[ > and ( |5.1 1 D give 


(5.11) 


/(Z0,£M rf £) < / [An - Rn-l](x 0 ,£)lp(d(;) < 1, 


where the second inequality is due to the definition of A n . Therefore, by Remark 4.19 
we get that xq, /-desires the centers in S 2 n B°. This gives an unstable pair for /, a 
contradiction. So suppose ip( S 2 D B°) = 0. We should have i/>(S 2 (T dBi) > 0. But ( [T9t 
gives />(<9 .Bi\S 2 ) = 0 and as before we get 

[ f{xo,£)v(d£) < I [An - R n -i](x 0 ,€)ip{dt;) < 1. 

JB i JSi 

So x 0 , f -desires the centers in S 2 n dBi, a contradiction. 

([TTJ) If x fully applies to £, i.e. A(x,£) = 1, then we have f s { x, £) = 1 — f?(x,£) by the 
definition of f s . Hence, the claim is a direct consequence of Q. 

din| Let A" be the set of sites x such that for if- a.e. £ the claim of © holds for (x, £). 
By © we have p(A c ) = 0. We prove that all sites in X satisfy the claim of ©f. Suppose 
£0 £ X and v Xo {f,t) > if Xo (f s ,t). Therefore, ip Xo (f s ,t) < 1. It follows that either x 0 is 
f s -unexhausted or f s (x 0 , •) is positive somewhere outside B t (x 0 ). Since f s is obtained 
by the site-optimal Gale-Shapley algorithm, in both cases x 0 has applied to all centers 


in B t (xq) with weight 1 (for the first case use Lemma 5.2 and for the second case 
note that x 0 has applied to some center outside the ball). The definition of X implies 
that f(x o,-) < fs(xo,-)> if-a.e. on B t (xo). Since v < if, we get Vx 0 (f,t) < if X0 (fa,t), a 
contradiction. 

©t By the right-continuity of p^(/, •) and •), it is enough to prove the claim for 

rational t. If this doesn't hold, we can find t £ [0, oo] and S C such that if(S) > 0 and 
p^(/, t) < t) for all £ £ 5. For arbitrary £ 0 £ S, since p > <p, the set of sites 


T «o : = i x e B t (£ 0 ) : f(x, £ 0 ) < / s (x,£ 0 )} 
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has positive p-measure. Moreover, we have p 5 ° (/, t) < 1 and /(•, £ 0 ) < 1 on Xj 0 . Therefore, 
Remark 4.19 gives that £ 0 , /-desires all points of We will show that £ 0 can be chosen 
such that some point of also /-desires £ 0 and contradiction follows. 

for p and ip, <[II} and Fubini's theorem on {(x, £) : £ G S, x £ T ,<} 


Assumption 


4.16 


(which is measurable) imply that there exists a site xo such that 

(a) ( |5.9P holds, 

(b) statement (|]]> is valid for x = x 0 and ip-a.e. £, 

(c) ip(C) > 0, where C = {/ £ S : x 0 £ T^j. 

In fact, these conditions are satisfied by p-a.e. xq. Consider the smallest ball B centered 
at xo (possibly with infinite radius) that contains C. Since f s {x 0 , •) > /(x 0 , •) > 0 on C, 
we get that x 0 has applied to the interior of B with full weight. Therefore, © implies 

f{xor) < fs{xo,-), f a.e. on B°. 

Case 1. ^(CnB°) > 0. By the above equation we get 

[ /OtnCM^) < [ fs(xo,Oi>{dO- 

J B" J B° 

Therefore f go /(x o,£)v(d£) < 1. So x 0 /-desires all centers in C fl B° by Remark 
which gives an unstable pair for /, a contradiction. 

Case 2. ip(C n B°) = 0. Therefore ip(C n dB) > 0. By fa} we get ip{8B\C) = 0. Since 
f s (x o, •) > /(xo, •) on C we get as before 

[ f(x Q ,Qip{d£) < [ f s (x 0 ,Qip(d£). 

J B J B 

Therefore v Xo (/, s) < 1, where s is the radius of B. So xq /-desires the centers in C fl dB 
by Remark |4.19| a contradiction again. 

□ 

Proof of Proposition \4.23\ Let / be an arbitrary stable constrained density for y and ip. 
Corollary |4.22| implies that for y-a.e. site x we have 


4.19 


1px{fs,t) =ip x {f,t), Vt £ [0, oo]. 


(5.12) 


We can use this equation for t < a s (x) (where a 3 (x) is the application radius of x in the 
site-optimal Gale-Shapley algorithm) and part (|n]l of Theorem |4.21| to obtain that for 
(y 8> ip)-a.e. (x, £), if \x — £| < a s (x), then / s (x, £) = f(x, £). Assumption 4.16 and ( |5. 1 2| ) 
for t = a s (x) imply that this is also valid for |x — £| = a s (x); i.e. for (y <8> ip)- a.e. (x,£), if 
\x — £| < a s (x), then f s (x, £) = /(x, £). Since f s (x, •) = 0 outside f? as (x) (x), we use ( |5. 1 2| ) 
for t = oo to obtain that (y <g> ip)- a.e. we have f s = /. This proves the claim. □ 

Proof of Theorem \4.24\ We should prove that almost surely any two stable constrained 
densities are equal except on a set with zero (8) 4^)-measure. Let f s and f c be 
the site-optimal and the center-optimal densities respectively. We take expectations in 
Corollary |4.22| and apply Lemma [5~4] to get 

A$E$ [4> 0 (fs, *)] > [Wo(/ c , t)] = A*E* [$°(/ c , t)] 

> A*E* [$°(/ 8 ,t)] = A*E$ [$o(/»,*)] • 

Therefore, all inequalities are indeed equality. Hence 


e $ [Mfs,t)-Mfct)} = o. 


Page 20j25] 






















Stable transports between stationary random measures 


By an argument similar to ( |5.3| ) we get 


E 



M/s,*) 




= 0 . 


Since the integrand is non-negative by Corollary |4.22| for a given t, we almost surely 
have for ^-almost every x G R d 


M/s,*) = M/c,*)- 

Considering this for rational t and using right-continuity of 4 ' x (f s , •) and 'b x (f c , •) we get 
that almost surely for 4> w -a.e. x we have 


’M/s,*) = M/c*), VfG[0,oo]. 


Now, for a sample w G ft such that the above equation holds for 4> w -a.e. x and Assump¬ 
tion [4716] holds for 4> w and 4/^, the claim is a direct consequence of Proposition |4.23| □ 


Proof of Lemma [4.'29l Remark 14. 2 81 easily implies that the Voronoi territory of £ is star¬ 


shaped with center £. As an example for non-convex territories, let i/> be half of the 
counting measure on the vertices of an equilateral triangle in the plane. Also, if ip = £+<5 0 , 
we see that the Voronoi territory of center 0 is not closed. 

Now, suppose the support of i/> is a discrete set {£ 1 , 62 , • • •}• Let D, be the Voronoi 
territory of a center £, ; . In this case. Remark 4.28| gives that x G I), if and only if 
| (x)) < 1. This easily implies that Di is closed. With the notations of Def¬ 


inition 


4.26 


let A be the set of points x G R d such that dB s ^ (x) contains more 
than one atom of i/j. In each component of A c , the set of atoms in B s m (x) is fixed. 
Therefore, all Voronoi territories are bound by A. Moreover, A C n ijA it j, where 
A h j := {x : C dB s ( x ) (a;)}. Each A, n] is contained in a hyperplane. Therefore, it 

suffices to prove that for each compact set K, only finitely many of the A, j's hit K. Let 
S D K be a compact set such that ip{S) > 1 (in the case ip(R d ) = 1 the claim is trivial). 
Each open ball B°^ (x) for x G I\ doesn't contain S and therefore it is contained in a 
compact set T which depends only on I\ and S. Now, if A, j hits K, then £, and ^ lie in 
T. It follows that such pairs (i,j) are finite and we are done. D 


Proof of Proposition \4.30\ Suppose the Voronoi territory C of center £ 0 is unbounded. 
According to Lemma [4.29| C is star-shaped and hence there is a half-line l starting at £0 

gives ip(B° x _ 


4.28 


which completely lies in C. For x G l, Remark ■ 
union over x G l, we find < 1, where H is an open half-space orthogonal to l, a 

contradiction. 

Now, suppose 4' is a stationary non-zero random measure. One can obtain from 
Lemma [2TT] that every half-space has infinite 4'-measure a.s. and the claim follows. □ 


Lemma 5.6. Let A C R d , a G lR. rf and r > 0. Let C be the set of points x such that 
\x — a\ < |a; - a'\ + r for all a' G A. C is bounded if and only if B r (a) is contained in the 
interior ofConv (A), where Conv (A) stands for the convex hull of A. 


Proof. One has C is star-shaped with center a. Hence, it is unbounded if and only if it 
contains a half-line {a + tv : t G [0, 00 )} for some unit vector v G R d . This is equivalent to 
A C {x : v ■ (x — a) < r} for some unit vector v. Equivalently, B r (a) % Conv (A)°. □ 

Proof ofTheorem \4.25l Using the ergodic decomposition theorem, we can assume that 
P is an ergodic measure for the family ( 0 s ) sgR d. 

([I) We call a center bad if its territory is unbounded. Let S' be the set of bad centers. 
Suppose 4'(S') 0 with positive probability. By ergodicity, this happens almost surely. 
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(b) A\ — R\ 


(c) / 


Figure 1: The supports of three functions in the setting of Example |6.2| which map 
[0, a] x [0, a] to {0,1}. The first axis stands for sites and the second axis stands for centers. 
The dashed lines are the lines with distance \ to the sides. The small segments in the 
third figure have slope — 0 or — ft where </> = 1+ ^ . 


Since F is flow-adapted, the restriction 'L' of 4/ to 5' is an ergodic random measure 
with positive intensity. Let ft £ S' and let C be the Voronoi territory of ft with respect 
to iH/'. We claim that the territory of ft is a subset of C and contradiction follows 
directly by applying Proposition |4.30| for If the claim is not true, let xq £ C such that 

F(xo,t;o) > 0. By Remark 


center ft £ B° XQ _ io{ (x 0 ) 


4.28 

such 


we get 'i 1 ' (x 0 )) > 2 > 1. Therefore, there is a bad 
that F( a ; 0 , ft) < 1. Now x 0 desires ft since it prefers ft 


over ft. Also ft desires x 0 since its territory is unbounded. So (a; 0 ,ft) is an unstable pair, 
a contradiction. 

|u) By stationarity, it is enough to show that for a deterministic point z, the union of 
the territories of the sites in B\ (z) has finite ^'-measure almost surely. Call a lattice point 
z £ 1 d bad if this property doesn’t hold for z. Suppose by contradiction that 0 is bad 
with positive probability. Then the set of bad lattice points form an ergodic simple point 
process on with positive intensity. Let Z be the set of bad lattice points, which contains 
zq := 0. Since every half-space contains some element of Z, we have Conv(Z) = Eft 
Therefore, we can find zi,...,z n G Z\{0} such that B 2 (zq) C Conv (z ±,..., z„). Let Co be 


the region defined in Lemma 5.6 for A = {zi,..., z n }, a = zq and r = 2, which is bounded 
by the claim of the lemma. For i = 1,... ,n let (ft be the set of points x £ Co that are 
closer to z,, than other points of A. Let T be the union of the territories of the sites in 
Bi (z 0 ). Since z 0 is a bad lattice point, we have 'T(X’) = oo. Since 'L(Co) < oo, there is 
i > 0 and a bounded Borel set D C T n Ci such that ’f'(.D) > 1. Since Zi is a bad lattice 
point, there is a site Xi £ B\ (z,) such that its territory contains some centers further 
away than all of the centers in D. Since 'L(D) > 1, we can find a center ft £ D such that 
F(xi,C) < 1. We claim that (a;*, ft) is an unstable pair. 

Since ft £ T, there is a site x 0 £ Bi (z 0 ) such that F(x 0 , ft) > 0. The fact that ft ^ C 0 
gives |ft — z 0 \ > |ft — z,| + 2. Therefore |ft - x 0 | > |ft — Xi\ and thus ft desires Xi. On the 
other hand, Xi desires ft since its territory contains centers further away than ft by the 
definition of x,. So (a:*, ft) is an unstable pair, a contradiction. □ 

6 Examples 

In the following examples we let Z d be the counting measure on X d ; i.e. 

^ = E 
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Figure 2: The territories of two sites x\ and X 2 and a center £1 in Example 6.4 The 


vertical solid lines represent Z x It. The dashed circles show the application radii of the 
two sites and the vertical bold lines show their territories. 


Although Zd is not a stationary measure on R d , in examples |6.4| and 6.5 we can translate 
the measures by a random uniform element of [0, l] d to obtain ergodic stationary random 
measures. Therefore, the claims of Theorem |4.14[ Theorem |4.8| and Theorem |4.24| are 
valid in these examples. 


Example 6.1. When ip = Cd and ip is a counting measure, the transport kernel given by 
the site-optimal density via Remark |4.3| coincides a.e. with the allocation presented in (3J. 
In this setting, definitions |4. 2 1 and |4. 12 1 also generalize the definitions of allocations and 
stability in Q. 


Example 6.2. Let ip and ip be both the Lebesgue measure restricted on [0, a], where 


a > |. 


Figure 1(a) and 1(b) illustrate the first stage of Algorithm 


4.4 


It can be seen that 


the site-optimal and the center-optimal densities agree a.e. with the function depicted 
in Figure [T(c)[ Therefore, we can use Proposition |4. 23 1 to see that this is the unique 
stable constrained density for p and ip. Note that the sites (resp. centers) that have 
distance less than ^(i _ i) to the boundary are unexhausted (resp. unsated), where <p is 
the golden number; i.e. it satisfies (p — \ = 1. 


Example 6.3. Let p 

and (ei one has 


C\ and ip be the counting measure on a subset A CZ. For x £ R 


fs(x,£) 


1 , I* -£| < \ 

0, otherwise 


It is easy to verify that f s is a balancing stable constrained density only when A = Z. 
Example 6.4. Let p = £ 2 and ip = Z\ ® C\. Define 


f 1, max{|a:i — £i|, \x2 — £2!) < 5, 

\ 0, otherwise. 


This function is a balancing constrained density, but is not stable since the center £ 0 = 


gives f s in one step, which is a balancing stable constrained density. For x = (xi,X2) and 
£ = (£i,£ 2 )> the site-optimal density is 

1, \X 2 ~ 6| < min{|, | — 2|xi — £i|}, 

0, otherwise. 

The territory of each center is a hexagon as illustrated in Figure [2] 



(0,0) and the site x 0 = (a, 0) desire each other for t < a < |. However, Algorithm 


4.4 
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Example 6.5. Let ip = 2 C\ and ip = C\ + Z\. For 0 < x < the site-optimal density is 


fs(x,0 


1, 0 < £ < 2x, 

1 — 2x, £ = 0, 

0, otherwise 


and for — \ < x < 0 we have 


r 1, 2x < £ < 0, 

fs(x,0=< 1-21*1, £ = 0, 

[ 0, otherwise. 


A similar equation holds for other values of r. By applying a random translation. 
Theorem |4.8|holds and f s is a balancing stable constrained density. 


Example 6.6. Let $ and be jointly stationary ergodic counting measures (i.e. simple 
point processes) in R, with positive and finite intensities A 4, and A4,. The random 
measures <i>' := ^-<I> and 4/' := ^-4/ have unit intensity, but there is no ($', 4d)-balancing 
allocation provided that 4 Z (note that ^ should be the number of pre-images of a 
center), ldowever. Theorem |4.8| shows that there is a flow-adapted balancing transport 
kernel between them. 


Example 6.7. In the setting of Example 6.6 let 4>" := 4>' x C d _ 1 and 4/" := 4>' x C d ~ 1 , 
which are stationary ergodic random measures on R d with unit intensity. Theorem 4.8 


shows the existence of a flow-adapted ($", ^"j-balancing transport kernel, but we claim 
that there is no such allocation. 

Let r be a flow-adapted (4)", 4’")-balancing allocation. By translation invariance, we 
get that for a site x, the vector t u (x) — x depends only on ni(x) and w, where 7ri is 
projection on the first coordinate. It follows easily that the allocation on R defined by 
x 7ri(r((*, 0))) is ($', H/'j-balancing, which is a contradiction by Example |6.6| 


Example 6.8. Let tp be the Lebesgue measure on (0, 00 ) and ip be the Lebesgue measure 
on (— 00 ,0). It is easy to see that f s (x, £) G {0,1} and f s {x, £) = 1 if and only if 
where [a] is the smallest integer not less than a. 
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